#include"main.h"
#include"ground_data.h"
#include"pipe_data.h"
#include"stationary_solver.h"

void main() {
  float depth = 10.0f;
  float width = 5.0f;
  int cell_q_y = 4;
  float min_size_x = 0.1f;
  float T_0 = 288.0f; // -273,3 = C
  float T_0_surface = 293.0f;
  float T_0_earth = 287.0f;
  float T_in = 303.0f;
  float T_out = 290.0f;
  GroundData* GroundData_ = new GroundData(cell_q_y, min_size_x, depth, width);
  float pipe_delta_x_inside = 0.2f;
  PipeData* PipeData_ = new PipeData(cell_q_y, depth, pipe_delta_x_inside);
  StationarySolver* StationarySolver_ = new StationarySolver(GroundData_, PipeData_);
  StationarySolver_->SetGroundInitialConditions(T_0, T_0_surface, T_0_earth);
  StationarySolver_->SetPipeInitialConditions(T_in, T_out);
  return;
}